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We show that the Josephson effect between triplet superconductors is sensitive to the relative 
orientation of the d- vectors across the junction. In addition, we point out that the temperature and 
angular dependence of the Josephson effect can help distinguish between different order parameter 
symmetries. We discuss coherent and incoherent tunneling processes, and the role of surface bound 
states within the framework of non-local lattice Bogoliubov-deGennes equations. 



Recent upper critical field and NMR experiments in (TMTST)2X [1-3], where X = CIO4, PFg have indicated that 
they may be unconventional triplet superconductors, while low temperature thermal conductivity experiments [4] 
suggest that the superconducting state is gapped. These experiments were inspired by early suggestions of triplet 
superconductivity in Bechgaard salts [5-8], and served as inspirations for the intense theoretical efforts [9-14] that 
followed. Lebed, Machida, and Ozaki (LMO) [11] proposed a "p-wave" triplet order parameter for (TMTSF)2PF6, 
where the d-vector had a strong component along the a direction, thus producing a strongly anisotropic spin sus- 
ceptibility with Xa <C xn and xv ~ Xn- A fully gapped singlet "d-wave" order parameter for (TMTSF)2C104 was 
proposed by Shimahara [12], while gapless triplet "f-wave" superconductivity for (TMTSF)2PF6 was proposed by 
Kuroki, Arita, and Aoki (KAA) [13]. Duncan, Vaccarella and Sa de Melo (DVS) [14] performed a detailed group the- 
oretical analysis and suggested that a weak spin-orbit fully gapped triplet "p^j-wave" order parameter, where Xa ~ Xn 
and Xb' ~ XN [10,15], would be a good candidate for superconductivity in Bechgaard salts. 

Here, we propose a test for the order parameter symmetry of Bechgaard salts which is based on the Josephson 
effect between two triplet superconductors (TS) separated by an insulating (I) barrier. Our suggestion differs from 
previous proposals invoking tests of the symmetry of the order parameter in the context of triplet heavy fermion 
superconductors [16,17], which relied on the Josephson effect between a singlet and a triplet superconductor. We 
show that the Josephson effect between two TS exists in the simplest case where the tunneling matrix element is 
spin-conserving (with or without time-reversal invariance), and that it depends on the relative orientation of the 
d-vectors of the TS. In addition, we show that the temperature and angular dependence of the Josephson effect can 
also help distinguish between different triplet states. We include coherent and incoherent tunneling processes and the 
effects of surface bound states (within the framework of non-local lattice Bogoliubov-deGennes equations). 

In anticipation of the existence of surface bound states in particular geometries of unconventional superconduc- 
tors [18,19], we prefer to use a real space representation of the left [L) and right (R) superconductors Hamiltonian 



i"j)c]_„^ (rj)cj- a,- (r'j) + 



E 



H.C. 



(1) 



where indices aj and f3j (with j = L,R) are spin labels, rj,r'j are position labels. Repeated greek indices indicate 



summation. The first term of Hj contains Tl{rj,r'j) 



) — fiSrj.r'j], where tj{rj,r'j) are transfer integrals 



confined to nearest neighbors and fj,j are the chemical potential. The second term of Hj contains the order parameter 



^i,a,-,ft-,5,,7,(i"j,r'j)(cj,5,(r'j)cj,7,(rj)), where y,-„^.,^^.,5,,7,(rj,r'j) is the two-body pairing 



matrix A 

interaction. The tunneling Hamiltonian connecting L and R superconductors is 



(2) 



where the matrix element for tunneling is Tq^q.^ (r/,, r/j) = Ts{rL,rii)6a 



- T„ (r^ , Tfl) • CTc^ ^OLR, with the first term 



Ts corresponding to spin-conserving tunnelling, and the second term T„ corresponding to spin-dependent tunnelling. 
The pair current through the junction is (to second order in the tunneling matrix element) 



n 



exp{iujjt) Ta^anii^L,rR)Ti3^^i3„ (r'i, r' r)P^^^^ (ri, r^, t'l, r' r, iuj,,) 



^' L^' R. 



(3) 



where iuon corresponds to Matsubara frequencies continued to (—eV/Ti + iS), and ujj = 2eV/h with V being the 



voltage applied across the junction. The tensor P^^^^ can be written in terms of the anomalous Green's matrices as 
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Pa!;fj^ii^L,rii,r'L,r'R,iuJn) = TJ2tu^ -F^^^^^^ (r^, r'^, ij/„)FK,„^0^(rfl, r';?, icj„ - iv^), where F,-,„^,3, (r^, r'^, iz/„) = 
T J dr exp{wmT)Fj^aji3j{'rj,r' j,t), with Fj,a^p^(rj,r' j,t),^ {Tcj^aj{rj,T)cj,p.{r'j,0)}. The general form of the order 
parameter matrix Aj^a.^f}.{rj ,r' j) — iAsij{rj,r'j) [(J2]a + (r^- , r'j ) • [crcr2]c /j. > where the first term corresponds 
to the singlet (pseudo-singlet) state and the second to the triplet (pseudo-triplet) state in the case of weak spin-orbit 
coupling (in the case of strong spin-orbit coupling). The order parameter matrix A^.a^.^^g^. (r^, r'j) satisfies the Pauli 
exclusion principle, since the function Asij(rj,r'j) which represents singlet pairing is symmetric under the exchange 
Yj ^ r'j. and the vector dj(rj,r'j) which represent triplet pairing is antisymmetric under the exchange Vj <-^- r'j. 
The reduced Hamiltonians on both sides of the junction can be separately diagonalized via the lattice Bogoliubov-de 

Gennes transformation Cj,a^(rj) = J2nj UNjaj{Yj)lNj '^*Njaj{^o WNj ' where UNja^ and VN^aj are two component 
spinors for each value of aj and satisfy the corresponding non-local Bogoliubov-deGennes (BdG) equation 

EN^UN^a^ {vj) = n{rj,r'j)uN,aj {r'j) + ^j,aj,i3, {rj,r'j)vNji3, {r'j) (4) 



-ENjVN^aj (rj) = J2 ""'jO^^iv^ai (r'j) + (""j, r'j)ujv^/3^ (r'j)- (5) 

This set of equations must be solved self-consistently together with the order parameter equation 

^j,a,-,/3,(rj-,r'j) = Jv,>,,ft,5,,7,(rj,r'j)X^(l-2/ArJ [y*N,^,{rj)uN,s,{r'j)-vl!.s^{rj)uN^^,{r'j)'\ , (6) 

where the two-body potential is written under the assumption of weak spin-orbit coupling as Vj^aj,i3j,5j,fj (i"j, r'j) = 
'^h,j{rj,r'j)Socj^^Sf}^^s. +2hj{rj,r'j)aoc^^j • 

From now on we will assume that the superconductors on either side of the junction are in (a) a triplet unitary state 
(time reversal invariant state in the bulk); (b) characterized by weak spin-orbit coupling; (c) the tunneling matrix 
element T'q.^^j^ {rj,r'j) is spin conserving, i. e., the tunnel barrier preserves spin, and thus it is not magnetically 
active; and (d) the dj-vector is weakly locked to the Cj axis by the weak spin-orbit coupling, a choice that is 
consistent with Knight shift experiments of Lee et. al. [3]. All these assumptions seem to be applicable to the 
Bechgaard salts. We model these quasi-one-dimcnsional salts via the bulk dispersion ej(kj) = —\tx.j\cos{kxjaj) — 
\tyj\ cos{kyjbj) — \tz,j\ cos{kz,jCj), where \tx,j\ ^ \ty,j\ ^ \iz,j\, corresponding to an orthorhombic crystal (I?2h group) 
with lattice constants a, b and c along the a, b and c axis respectively. In the D2h point group all representations 
are one dimensional and non-degenerate [14], which means that the d- vector in momentum space for unitary triplet 
states in the weak spin-orbit coupling limit is characterized by one of the four states: (1) ^Aiu{a), with dj(k) = 
fijAf^^^^^jXjYjZj {'\Uy," state); (2) •'Si„(a), with dj(k) = ^Ap^jZ^- ("p/' state); (3) ^B2u{a), with dj(k) = 
fjjAp^^jYj {''py" state); (4) '''i?3,j(a), with dj(k) = fjjAp^,,jXj {"p,,." state). Since, the Fermi surface touches the 
Brillouin zone boundaries the functions Xj, Yj, and Zj need to be periodic and can be chosen to be Xj = sin {k^jaj), 
Yj = sin {kyjbj), and Zj = sin {kz.jCj). The unit vector fjj defines the direction of dj(k). Next, we discuss the 
Josephson effect in two situations, first neglecting, and second including surface bound state effects. 

The Josephson effect neglecting surface bound states: In this approximation only the contribution from scattering 
states with well defined momentum is included. Under this assumption the BdG amplitudes become 'UjVjaj (rj) = 
exp{i]s.j ■ rj)ui^. / W ; VNjOj {rj) = exp(ikj ■ r j)v]s.j / W , where quantum indices Nj are represented by momentum 
kj and discrete index rij = lj,2j. In this case, the expression for Js{V,T) defined in Eq. (3) becomes 



WT) = | 



cos((/)) Y Im{W{kL,kR,iuJn)} + sin{(j)) Y Re{W{kL,kR,iujn)} 



(7) 



after summations over spin indices oij-,(ij, indices rij, Matsubara frequencies ivm, and the global U(l) gauge 
transformation dj djexp(z(/)j) arc performed. Here, <f) = <j)R — 4'l + ^jt, and W(kL, k/f, ia;„) = 

2dl{kL)-dR{kR)Q{kL,kR)P{kL,kR,iujn), with Q{kL,kR) = r,(kL, kH)r,(-kL, -k^), and ia;„ ^ -eV/U + iS. No- 
tice that the dot product d|^(ki:,) • dii(kii) appears in Js{y,T), and thus Jsiy,T) is very sensitive to the relative 
orientation between d|^(ki) and di^(ki^) [20]. For instance, if d^(ki) ± d7j(k/j) then Js{V,T) vanishes identi- 
cally for any V and T. Furthermore Js{V,T) (for fixed V and T) changes sign depending if the vectors dj^(ki) 
and di{(ki{) are aligned or anti-aligned. It is this sensitivity to the relative orientation of the d- vectors across 
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the junction that makes the Josephson effect between TS a crucial test for triplet superconductivity itself. The 

function P(kL,kR,lLOn) = TJ2i,^ [i^^m? + [i^^' 

^i^m = (2m + l)TTkBT. The excitation energies are i?kj 
identically if VF(kL, k_R, zw„ —eV/h + i6) is odd under inversion for Icl 



i?k is even under inversion k, 



kj, and P(kL, kij, iLu„ 



^VmY + -E'kn] ^ ^^'^ Matsubara energies 

p. Notice that Js(y,r) also vanishes 

kL or k_R, ^ — kij. Since dj(kj) is odd, 
-eV/fi + is even under k^ — k/, or k/j — kij, 



e + |d,(k,) 



a non-vanishing JgiV^T) requires Q(kL,kij) to have contributions which are odd in both k^ and kj^. Therefore, 
Js(y,r) depends crucially on the magnitude and symmetry properties of (5(ki,kjj) defined above, which contains 
the tunneling matrix element T's(ki,kjj). Since, all the Z?2/i point group representations are one dimensional and 
non-degenerate [14], can be expanded as Ts(ki, kj^) = "^^TlTr '^TlTr(}^l,^r)'4'^^{^l)'^^^0^r)^ in the case of weak 
spin-orbit coupling. Here il}^^{k.j) are the basis functions of of the D2h (orthorhombic) point group [14]. In the 
particular case of time reversal symmetry rs(— k^,— kj^) ~ r*(ki,k7j) and QikL,^.^) — |Ts(ki, ki^)!^ > 0. 



1 

0.8 

© 

0.6 

H 

X 0.4 
0.2 














1 






0.5 













■0.5 






■1 




\ \ ■ 


45 90 135 180 ''•■^ 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

T/T, 

FIG. 1. Plot of Js{T)/Js{0) versus T/Tc- Inset shows Js{0,aLR)/ Js{0,0) versus q_l_r at T = 0. Several a axis tunneling 
processes are illustrated for "pa;" and "py" symmetries. Incoherent processes: "pa;" with (a) Q'*"'^-' oc XlXr (solid line) and 
(b) Q^*""^ oc lilfl (trivially zero); with (c) Q^'""^ oc YlYr (dashed line) and (d) Q'^^^ oc 1^1^ (trivially zero). Coherent 
processes: "p^:" with (a) oc XlXr (dotted line) and (b) Q'^""''^ oc IlIr (trivially zero); with (c) Q'-""''^ oc YlYr 

(dot-dashed line) and (d) Q^'^°'^^ oc IlIj; (double-dotted line). The L and R superconductors are assumed to be identical 
quarter-filled systems, with Tc — 1.5 K, and parameters t^ = 5800 K, ty — 1226 K, t^ — 48K. 



For the weak spin-orbit "p^", "p^"' "Pz" ^^fxyz" states, where the direction of the d-vectors is independent of 
momentum, a simple expression for Js (T) at zero bias can be obtained 



jW{T) = J(^)(r) X {fiL ■ m) X sin(0), 



(8) 



where j(^)(r) = (2e/?i)Ai(T)Ai^(r)S'j^-^2(T), (l)^4>j^~ 0^, and the label A identifies incoherent (A = inc) or coher- 
ent (A = coh) processes. Here, sfl{T) = i^^^g(^)(ki, kfl,)V'^^>"(kL)V'^«^"(k,j)P(ki, k.^), where P(ki,kfl,) = 

4TE™=o(?i'^™ + ElJ-\h^vl + El^)-\ and ^'^■"(kL) ^ X, Y, Z or XYZ. Notice that t}l ■ m = cos(0Lij), 
where cos{9lr) is the angle between the two d-vectors, and describes a polarization angle just like the Malus's law 
of electric polarization. Js{T) depends crucially on the tunneling matrix elements (5(ki,kjj). For purely incoherent 

processes the only terms that contribute to ji™'''' (T) come from Q^*"'^) = 2T}^^'^T^^^" - 2T}^^'^T^^^'^ , where fj are 

appropriate odd representations of the dj-vectors. For the "p^^" symmetry V'^^'"(kL) = Xl, and ip^"''^{'kR) = Xr, 



thus (a) Qp"p}{'kL,'i<-R) = Qp^p^XlXr produces a non-vanishing Js{T), while (b) Qpjp j (k^ , k^J = Qp^p^Il^r 
produces a trivially vanishing Js{T). Similarly for the "pj^" symmetry -(/''"^'"(k^) = Yl, and ■(/''""'"(ki?) = 1^, 



(inc) / 
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thus (c) Q^p^p^ikL^kR) = Qp^pYlYr produces a non- vanishing Js{T), while (d) Q^"p';] (Icl, k^) = Qp^^p^lL^R pro- 
duces a trivially vanishing Js{T). For comparison purposes, we analyse next four coherent tunneling processes 
allowed by group theory corresponding to matrix elements (a) Q^°p] {k.L,'kR) = Qp^p^-^L-^ij^ky^ .ky^ , and (b) 
Qp°pJO^L,'k.R) = Qp^p^li,lij5k||^,k||^, for the V" symmetry; (c) Qp^pJ ik.L,'k-R) = Qp^p^ YlXr, (^kn^ .k,,^ , and (d) 
Qp!^pJO^l,^r) = QpaPylLlfli^kii^.kii^, for the "p,^" symmetry, where the momentum parallel to the junction is con- 
served. In Fig. 1 we show both the temperature and angular dependence of ji^"* for the case of a axis tunneling, 
assummg that dj II Cj (j = L,R). /s^\t) / Ji^\0) is shown for Q!L = 0, where all L and R axes coincide. In the inset, 
we show Js^\o, ulr)/ ji^\o,0) for T = 0. For coherent processes ky^ = ky^ and the parallel momenta are related 
by kyff = cos{aLR)ky^ — sm{aLR)kz^, and fc^^ = sm{aLR)ky^ + cos{aLR)kz^, where ulr is the angle between k||j. 
and k||^ . For a axis tunneling Olr = a^R. Notice in the inset of Fig. 1 that the Josephson current does not change 
sign for the "pj^" with (5*^^°'*' oc IlIh (double-dotted line). This is a direct manifestation of the polarization effect of 
the d vector. It is important to emphasize that both the temperature and the angular dependencies of Ja(T, a^ij) 
can help distinguish different symmetries, and different (coherent or incoherent) tunnelling processes (see Fig. 1), as 
seen experimentally for c axis tunneling of cuprate superconductors [21]. The approach used for a axis tunneling is 
good only in the cases of "p^" or "px" symmetry but not for "p^" or "Pxyz'N since the d-vector changes sign upon 
reflection at the interface, and leads to low energy surface bound states to be discussed next. 

The Josephson effect including surface hound states: To obtain Js{T) for the "p^" symmetry in the case of a 
axis tunneling it is crucial to solve Eqs. (4), (5) and (6) self-consistently. We solve the non-local lattice BdG 
equations and compare our results with standard local quasiclassical continuum approximations [22-24]. First, we take 
advantage of the translational invariance along the direction parallel to the interface and write the BdG amplitudes as 
UNjaji^j) = exp(ikj|| • rj\\)\^^^ jf,^^^{rj_)/VV ; VN^aji^^j) = exp(ikj|| • rj\\)v^^^^ j^^^^,{rj_)/VV, where quantum indices 
Nj are represented by momentum k^ and discrete index nj — Ij , 2j . After a Fourier transformation in the parallel 
coordinates the lattice BdG equations (4) and (5) become one-dimensional in the perpendicular coordinates rj_,r'j_. 
We choose hard boundary conditions, where the BdG amplitudes vanish when ry^ = ry^ = (at the center of the 

insulating barrier). Second, only the triplet channel component T^^*^ (r^ , r^- ) of the general weak spin-orbit coupling 

interaction l^',c«j,(3j,5j,7j (rj, r'j) is considered. The triplet component is Vj*\rj,r'j) = Iij{rj,r'j) + l2j{rj,r'j), and 

it is assumed that V^*-*-* (r j , ) has an on-site (r^ = r^) repulsion Uj and nearest neighbor attractions {r'j = Yj + a,;) 
Vj^xi Vj^y, and Vj^z only. In addition, we assume that the direction of the d(r,r') is independent of position and 
points along the c axis. Furthermore, we choose the spin quantization axis to be along the direction of d(r, r'). 
Under these conditions the order parameter matrix Aj^Q,^^^^. (r^ , r'j) for the "p^;" symmetry vanishes identically when 
r'j = Tj, and has only off-diagonal matrix elements Aj,1-^,j^^.(rj,r'j) = Vj^^ [(c^.t, (r'j)ci,j,, (fj)) + {cj,ij{i^' j)cj,ij{rj))] /2 
and Aj^i.^^.{rj,r'j) = Aj^-^.^i.{rj,r'j) when r'j = rj + aja. Js{T) is calculated numerically using Eq. (3) with 

TaLaniJ^LjTCR) = Ts{rL,rR)SaLa,i- For definiteness and simpHcity we also assume that Ts(rL,r_R) ~ TSr^^^^r^^^ for 
rL± = —d/2 and rR± = +d/2; Ts{rL,'rR) = otherwise. Here, d is the separation between the last layer of the L 
and the first layer of the R superconductor, and corresponds to the thickness of the insulating layer. This simple 
choice guarantees that the parallel momentum ky is conserved across the junction, and thus correspond to a coherent 
process. Since the tunneling matrix elements connect mostly states that have appreciably large BdG amplitudes near 
the junction, we use a mixed representation of Eq. (3) expressed in terms of the parallel momenta k||^ = '^^IIk = '^^ll 
and the perpendicular coordinates and to calculate numerically Js{T). In Fig. 2 we compare our results for 
Js(r) using the non-local lattice BdG equations and the local quasiclassical approaches [22-24]. 
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FIG. 2. Plots of Js{T)/ Jrci{0) versus T/Tc for for coherent "p^" ("pz") a (c) axis tunneling, where Jrcf(O) is the Landau 
critical current for a singlet s-wave superconductor with the same Tc- Quasiclassical results are indicated by dashed (dotted) 
lines for "px" ("pz"), while the non-local BdG results are represented by solid (double-dotted) lines for "px" ("pz")i for the 
parameters of Fig 1. 



The local quasiclassical approximation is only strictly valid when fcFi^_L ^ oo, but it leads to zero energy bound 
states in the case where the quasiclassical order parameter A(k, r) changes sign upon r — r [22-24] (as in the case 
of the a axis tunneling for the "pa:" symmetry). However, in the non-local lattice BdG equations (4) and (5) there are 
only finite energy bound states. The non-existence of zero energy bound states in the non-local lattice BdG equations 
can be understood as follows. Because the order parameter is non-local, it can be described in terms of the center 
of mass R = [r -I- r']/2 and relative r^ei = [r' — r] coordinates. Near the surface these two coordinates are entangled 
and thus lift the zero energy bound states degeneracy found in the local theory. This produces, perturbatively, finite 
energy bound states of the order |Ao|/7_l where ^± = f±/a_L, with and a± being the coherence and unit cell 
lengths, respectively. The finiteness of the energy of the bound states cuts off the low temperature 1/T divergence 
of Js{T) calculated in the quasiclassical approach at small transparencies [22-24]. For a axis tunneling and the "px" 
symmetry this amounts to a small correction to the quasiclassical results, as "fx = ^x/dx ~ 106 for the Bechgaard 
salts and Ap^p^ — 3.73 K, leading to the lowest bound state energy to be approximately T*^ = 35.2 mK. Which 
means that for T < T*^ the quasiclassical results are not reliable (see Fig. 2). The quasiclassical results fail in a more 
dramatic way, when we consider c axis tunneling for the "pz" symmetry. In this case £,z ~ 21.5 A, — 13.5 A, 
Iz = Cz/oz ~ 1-59 and Ag^p^ = 3.21 K, resulting in T*^ = 2.02 K, which is larger than the critical temperature 
Tc — 1.5 K used here. Thus, perturbative corrections to quasiclassical zero energy bound states are enormous, and 
the quasiclassical approximation can not be used except for T r^T^ (see Fig. 2). 

In summary, we have discussed the Josephson current J, in TS/I/TS junctions, and indicated that the triplet 
nature of the order parameter for Bechgaard salts could be tested in such an experiment. We emphasized that the 
temperature and angular dependence of Jg can help distinguish between different triplet order parameter symmetries, 
and we discussed the role of surface bound states within the non-local lattice BdG equations. We would like to thank 
NSF (Grant No. DMR-980311) for financial support. 
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